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FRACTIONAL RELEASE OF A TRACER ELEMENT 
THROUGH A MOVING BOUNDARY 


Stephen D. Lowe” 


Lawrence Radiation Laboratory 
University of California 
Berkeley, California 


October 30, 1963 


ABSTRACT 


A full solution for the fractional release of a tracer element 
through a moving boundary has been obtained for the cases of a slab 
and a sphere. In both cases, the initial distribution of the tracer is 
assumed to be uniform throughout the body. For a constant rate of 
boundary motion, the full solution is valid for all times over the range 
from zero to complete evaporation of the body, and this solution may 
be applied to all finite rates of boundary motion. 

By applying certain limitations to the time range and to the total 
amount of boundary motion, the full solutions are reduced to approxi- 
mate forms that are presented graphically as a family of curves. 

For a more limited range of application, i. e., for small values 
of the quantity 5 (5) (=) 1/2 (<0.2), where b is the rate of boundary 
motion, t is the time, D is the diffusion constant, and a is a charac- 
teristic dimension of the body (half-width for the slab and radius for 
the sphere), simplified expressions are developed for the fractional 
release of the tracer, f, which may be used to more accurately deter- 


mine the diffusion coefficient, D. 


These E expressions are 


m C» 1/2 | 1 (2) lab, 
7% 2 \a 


and i 


1/2 , + 5 (>) for the sphere. 


Ze (25 
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Pe ERODUCTEON 
A. Basis for the Problem 


During his investigation of the diffusion of xenon in uranium 
monocarbide, Shaked noted the slope of his plot of the fractional re- 
lease, f vs (t) 1/2 increased noticeably during a high-temperature 
EE He postulated that evaporation of the uranium monocarbide 
was responsible for the increase. 

A similar increase in the slope of the f vs (1) 1/2 plot had been 
reported previously for experiments in the temperature range, 1600 
[02220076 S 

With the increasing interest in high-temperature applications 
for uranium ftuels, 1n thermionics, for example, itis considered of 
Eniieresttosınvestigate the release of a fission product from the fuel 
element due to the additional process of evaporation occurring simul- 


taneously with the process of diffusion. 


B. Discussion of the Problem 


Several investigators have reported diffusion constants for 
xenon-133 and krypton-85 in uranium dioxide." A lesser number 
have reported results for the diffusion constant of xenon-133 in 
uranium monocarbide. * In the ''conventional' analysis, one assumes 
AUO initial concentration of the tracer element in the geometri- 
cal body whose boundaries are fixed with time. Evaporation of the 
body is discussed in a general manner, but the effects are minimized 


experimentally. is 


C. Statement of the Problem 


This report is an investigation of the effect of evaporation upon 
the fractional release of a radioactive tracer element from a fuel 
body. As such, the decay of the tracer is neglected for the purposes 
of this report (i.e., A=0), although this effect can and should be 
taken into account for experiments involving long postirradiation 
anneal time. 

An expression for the combined fractional release due to both 


diffusion and evaporation (i.e., diffusion through a moving boundary) 


is obtained. This result is compared to the fractional release derived 
for the case of ''pure diffusion'' (i. e., diffusion through a stationary 
boundary). 

The analysis is performed for the case of the finite slab and the 
case of the finite sphere. Itis assumed that the initial concentration 
of the tracer element in the body is uniform and that the concentration 
of the tracer element at the boundary (moving or stationary) of the 


body vanishes. 


II. THE SLAB PROBLEM 
A. Concentration of the Tracer 


One seeks a solution for the diffusion equation 


ut) e — (a-bt) « x € a-bt (IL-4) 
> 
such that CX O= C (II-4a) 


and E l=(a- bt) tJ = O et < a/b, (II-4b) 


where b is the rate of movement of the boundary that may be obtained 
from kinetic theory if one assumes evaporation in a vacuum ofa 
material of known vapor pressure. S The assumption of a constant 
evaporation rate ignores the possible effect upon vapor pressure, and 
hence upon the evaporation rate, if there 1s a change of chemical com- 
position of the solid surface accompanying evaporation. At time 
t-a/b, the sample will have completely sublimed. 

Chambré has shown? that the solution to Eq. (II-1) that satisfies 


an initial even function f(£) and Eq. (II- 1b) may be expressed as 


C t) 7 og J FG) exe RE 


eer [fep [b G*79]5 t5 4] 


(II -2) 


where for the case under study 


f(s)2 C, (I -2a) 
and S(x,tis)= exp LD lot Es exp |- Cosy |. (IL-3) 


4 DC 


Cu 


Making the substitutions into Eq. (II-2), one obtains the result- 


ing expression for C(x, t): 


Ct)" zip [fen Feld 
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Performing the integration over & (see Appendix D), one may show 


(11-4) 
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Making the definitions 


e DE (IL- 7) 
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and 


_ ba 
E D ) (II -8) 


where 0o«& gt «i, 


and E. into Eq. (II-6), one obtains 


C T)= Co [et ( 32%) ret (2 5.1 
" Ci) exp) ne in(-89- a ferf An( I~ Le Ya 


-erf Zrll et) leva + $ (IL-9) 
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By using Eq. (C-4) and setting B=0 (i. e., no evaporation or boundary 








motion), Eq. (11-9) becomes 


Ch, T) = Co- £e [erre (e) erre 1 
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This may be rewritten, after combining terms, as 


Ca?) = C. Ss Sien [er lr Da + eng Garber | eet 


Equation (Il-11) may be recognized as one equivalent to Eq. (A-6), 
which shows that the expression for the tracer concentration does re- 
duce to that for the stationary boundary case when no evaporation is 


considered, 


E 


B. Total Fractional Release 


The total release of the tracer element due to the combined ef- 
fects of diffusion and evaporation may be obtained by integrating the 
concentration at any time (7) across the volume of the slab, dividing 
the result by the total initial amount of tracer present, and subtract- 


ing this quotient from å. Thus, 





f= Q- Q(T) 21-80 (SEN 
er SE 
where 
= 4 (2aC.) (II- 42a) 
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Performing the integration and combining, one has 
f= e + TR [terfe La - verte E | 
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"tU rate exp m lange} (ee pe! E 20 | 


ZU (IT-14) 
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GER 
For r=0(t=0), the fractional release, given by Eq. (II-14), is 
zero, as is to be expected; and for 87-1 (or complete evaporation), 


the expression becomes 


+ = 1 + TR herk z FE Solla Le rl ze, 
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| Se 
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which is equal to 1, after cancellation of equivalent terms, and cor- 
responds to the total release of the tracer. 


Setting B=0 (no evaporation), Eq. (II-14) becomes 


T= T C "terc tn 
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This may be reduced to 
/ d = : 
{= 3 [7 ROSA verte TUA J (IL-16) 
nz| 


which may be recognized as equivalent to Eq. (A-8), the fractional re- 
lease under conditions of a stationary boundary. 
C. Approximate Expression for the Fractional Release 


For the usual values of 7 encountered in fission gas-release 
experiments (7 < 0.04), the expression obtained for the fractional re- 


! 
lease in Eq. (I1 - 44) may be immediately reduced to 


f* gt * U^ verte £5, 
Mae ntm m-l 
ORC Sp Ergo] ener]. ios Creo! 


(IL-17) 
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If the further limitation is made that B7<0.5, Eq. (II-17) may be 


reduced further to 


la. 
f* gt t “verte £5, 
11-18) 


En m~] l 
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xs 


of which the n=1 term is of the greatest significance, as follows: 


Pa pr + TA iiec t. 
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By using the relations (C-13) through (C-18), Ea. (II -19) be - 


comes 
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(IL-20) 
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4/2 4/2 
For large BT (i. e., Pr — >2.5), Eq. (I-20) reduces to f=ßr, 


which indicates that the release of the tracer is due solely to evapora- 


tion. 


For 8-20, Eq. (IL-20) becomes 


i = TUA , (11-24) 
Ro 


the approximate expression for fractional release from a slab with 
stationary boundaries. 

Figures II.1 and II.2 show the fractional release from a slab vs 
[time] 1/2 (73/2) for various rates of evaporation (8). As a reference, 
the stationary -boundary plot is given, and in each figure it is the 
lower line. Figure II.1 is for the longer periods of time, whereas 
Fig. II.2 represents shorter periods of time, 

Figures II.3 and II.4 demonstrate the contribution of the individ- 
ual terms that make up the expression [Eq. (II-20)] for f. Figure 
II.3 represents an evaporation rate of 8 - 10 for times less than 
7=0.01, and Fig. II.4 covers the time span less than 7=0.0001 fora 


higher evaporation rate of 8-50. 
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Fig. Il.1. Fractional release vs y 1/2 for various rates 
oí evaporation. 
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Fig. II.3. | Components of fractional-release curve. 
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Fig. 11.4. Components of fractional-release curve. 
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Ill. THE SPHERE PROBEEM 
A. Concentration of the Tracer 


A solution for the diffusion equation is sought. Let 


P avi) DI ge £2 2c S NOS (III - 4) 
such that C Cr. 0) - (s (III -4a) 
and C(a-bt,t)=O for O4t <4, (to) 


where b, again, is the rate of boundary movement. At time t=a/b, 
the sample will have been completely sublimed. 
Equation (III[-1) may be reduced to that for the slab problem by 


making the substitution 


u(rt)= rc Aa), (ii 
subject to the conditions 

ROEE E (III-2a) 
SE U(a-Lt t)= O for DATA, (IIL -2b) 


Since C(r, t) is finite at interior points, itis necessary to define 
ul Sech = 0 (III -2c) 


m Ee has shown the solution of Eq. (III-1), which satisfies an 


initial even function f(Ẹ) and Eq. (III-1b), to be 
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where, for this problem 
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By making the substitutions into Eq. (III-3), the resulting expression 
r Cr, t) is 
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By performing the integration over & (see Appendix E), Eq. (III-5) 


may be shown to ^. 
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Using the definitions (II-7) and (I1-8), then substituting into (III -6), 
one obtains 
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Setting 8-0 (i.e., no evaporation or boundary movement) in 


Eq. (III-7), it follows that 
C (E t) = B —- Loa Terfe T. — erke ie 
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The ierfc terms reduce to zero, and after some recombination, 


Eq. (III-8) becomes 
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This may be recognized as equivalent to Eq. (B-6), which shows the 
expression for the tracer concentration to reduce to that for the 


stationary-boundary case when no evaporation is considered. 


B. Total Fractional Release 


The total release of the tracer element due to the combined ef- 
fects of diffusion and evaporation may be obtained by integrating the 
concentration at any time (7) throughout the volume of the sphere, 
dividing this result by the total initial amount of tracer present, and 


Or e Ane this quotient from 1. Hence, 
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Integrating and ne the terms in (III-11), one finds 
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For r=0(t=0), the fractional release, as given by Eq. (III-12), is 


zero, as expected; and for 87-1 (or complete evaporation) 
pression becomes 
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which is equal to 1, after cancellation of equivalent terms, and cor- 
responds to total release of the tracer. 


Setting B - 0 (no evaporation), Eq. (III-42) becomes 
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This may be reduced to 
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S 
which may be recognized as equivalent to Eq. (B-8), the fractional re- 


lease under conditions of a stationary boundary. 


C. Approximate Expression for the Fractional Release 


For the usual values of T encountered in fission-gas-release 
experiments (7 < 0.01), the expression obtained for the fractional re- 


lease in Eq. (III-42) may be immediately reduced to 
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By making the further limitation that 67 < 0.5, one obtains 
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of which the n=1 term is of the greatest significance, as follows: 


f= 36t - Berl +(@t) + ST" (1-prierk t, 


GL L "erte Ze, + 6 rg Uert Zen, 
- I2 t i*erfe Et E 


+37 1-0) 261)” Get] certe EL 


+6 TEN mare)" l ( Morte a 


- 6 (I- PILLE I)” (ac ^)" " L T NE 
— 2 CE C m-2X2 pc ^Y imei Ze: | 


By using the relations (C-13) through (C-19), Eq. (III-47) becomes 
f= 3ec(i-evt) +t) - 3T - ar. 
+67" SR £t) eric BE" e (C) Cerfe -- 


(IIL-47) 


mod 


+I (e)z UTEE T erte BC" 
(III - 18) 


LG TU-SE etar P (Merle SE 


+ G ee (m C4 ) "LAL KS a 1 p 


oc )" EN oe 
+12 T Zi (muy) (e T SA BE. 


The last three summations will be small compared to the first terms: 


therefore, f will be approximated by 


£- Set (i-e) (6 - s c- 3c) £t (Ioan) lef 8E 
: zm - Sal + en (4) "(ace)" Verte ex^]. 


(III -19) 


For large B ile E > 2.5), Eq. (III -19) will reduce to 


f» Sgc(i-go) 6 - 8t - Stt (1-09), 


puscbuudicstes that the release of the tracer is due primarily to 
evaporation. 


mom = 0, Eg. (I-19) becomes 


VE CE 
ae T oR (IIL-20) 
the approximate expression for fractional release from a sphere whose 
boundaries are stationary. 
If a further restriction is made, such that ee 0701 (or 
tr <0.0001), then those terms in Eq. (III-19) having T as a coefficient 


(excepting BT terms) may be neglected, and one obtains 


2A = 


$e scm TET l, u (IIT-24) 
Schier) Rierk &- +2 (et^ inb ge^ T. 





1/2 
For large -- , Eq.(III-21) reduces to 


+= 37 (1-67) +(87)- 


which indicates the release of the tracer to be due solely to evapora- 
tion. 


For 8-20, Eq. (III-241) becomes 


NEN 773 
f- 2 C (IET 2293 


the more familiar expression for fractional release from a sphere 
with stationary boundaries. 

Figures III.4 and III.2 display the fractional release from a 
sphere vs [time] 1/2 (71/2) for various rates of evaporation (8). As 
a reference, the release for the case of a stationary boundary (f= 0) 
is given for each figure; see line (1). Figure III.1 (note the dotted 
line) in addition, shows the effect of neglecting the 37 term in 
Eq. (III-20) for longer periods of time. Figure III.2 represents the 


release of the tracer element for shorter time periods. 


ELSE 








AO O 1/2 
e 
(1) B= 0 (no evaporation) (3) B= 40 
. 6 4/2 : 
e 11/2 37 (4) B= 20 
Se = 20 
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O 0.02 0.04 0.06 0.08 0.10 
T 1/2 


MU-32929 


Fig. Il.i. Fractional release vs 71/2 for various rates 


of evaporation. 


225 - 


(1) 6 =  O(no evaporation) (5) B = 500 
(2) B= 20 (6) B = 1000 
(3) B = 100 (7) B = 2000 
(4) B = 200 
(7) (6) 0.06 
(5) 
(4) f 
(3) 0.04 
(2) 





0.02 
O 
O 0.002 0.004 0.006 0.008 0.01 
7 v2 
MU-32930 
; 112 
Fig. IIL.2.  Fractional release vs 7 for various rates 


of evaporation. 
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IV. DISCUSSION 
A. General Results 


For all values of 870, the fractional release with time is in- 
creased over the corresponding value for f when ß=0. For any value 
of time (or 7), f will increase, for both the slab and sphere cases, 
as B increases. With smaller values of B, say 8«20, the curve of 
1/2 1/2 < 9 04, 
which indicates that f is proportional (or very nearly so) to 71/2, 


1/2, 0.1, provideaw rn 


f vs T very closely approximates a straight line for T 


Ihe same statement may be made for 7 


b. Initial Release 


This theory does not account for the large early release rate in 


176 


postanneal experiments noted by several observers. Insehese 
cases, the slope of the f vs ¿1/2 curve decreased with increasing 


time until a region of more or less constant slope was attained. 


C. Spherical Case--Stationary Boundary 


Figure IIL.4 illustrates the effect upon f of neglecting the second 
term of Eq. (III-20). For values of u: 0.04, the difference in 
volwessor si resulting from lack of consideration of the second term on 
Eq. (III-20) is negligible. For larger values of De) the fractional 
pelease predicted by Eq. (III-22) becomes increasingly greater tham 
that predicted by Eq. (III-20), and differs by 0.03 (x 10%) at 71/2 - 0.1. 
Since both Eq. (III-20) and Eq. (III-22) appear to be plotted as straight 
lines, an error of 17% in the value of D results for the above value of 
7 from neglecting this term. It is noted that Eq. (III-22) closely 


represents the release curve resulting from boundary motion for 


which ß=2. 


D. Spherical Case--Moving Boundary 


For values of 71/2 < 0.04, Eq. (Hl=21) may be used in preferente 
toma. (111-49) for all values of Br. If Br (remember QT independent 
of D and is bounded between O and 1) is approx40.2, then the term 


(pr)? 


may also be neglected. 


20. 


For large values of gr V? the fractional release 1s approxi; 


1 
2 
mated by 3ßr(1-Pr) +(pr)3, 


This may occur at short times for very rapid boundary motion (8), or 


which is a cubic equation in time (or 7). 


for very long periods of time (7). Neither of these cases, however, 


is useful in making a determination of D. 


4 e 


For small values of 7 B , pertaining both to relatively short 


periods of time and small rates of boundary motion, the fractional re- 


1/2 


lease is directly proportional (or very nearly so) to T and there- 


1/2 


fore to 7 It should be possible, therefore, to develop a simple ex- 


pression that may be used to determine D with greater accuracy than 


871 2 B71/2 
would the use of Eq. (III-22). For small CHO "n < 0.2) ; hence, 
Eq. (III-241) may be approximated by 
A 
A 
E A be) E GEN, IV -4 
E A E ee (IV -4) 


With Eq. (IV-1), one can predict the release curve to be para- 
bolic in shape when f is plotted vs time to the one-half power; the 
curvature depends upon b. The curve consists of a straight-line com- 
ponent and a parabolic component, which is the evaporation correction. 

Figure IV.1 is the graphical presentation of Eq. (IV-1) for those 
values of f that do not differ by more than approximately 0.004 from 
the more accurate expression for f. 

By assuming that b (the rate of boundary movement) has been 
determined, a corrected value of the fractional release may be ob- 
tained by subtracting the quantity - (E) from the actual fractional 
release. This corrected value of f may then be used with Eq. (III-22) 


to determine D, 


E. Slab Case--Moving Boundary 


For large values of a the fractional release is approx1- 
mated by T--the percent of the body that has sublimed. 

For small values of a , vo 0.2) , a Simplified expres- 
sion for the fractional release may be developed in a manner similar 


to that for the spherical case. This is written as 





29% 


{= > BT + OR 73/2 
tu? 
(1) p= O (no evaporation) (4) B= 200 
wann = 20 (5) B= 500 
(3) p = 100 (6) B - 1000 
0.06 
f 
0.04 
(ap^ C» 
(5) 0.02 
(6) A 
O 
d 0.002 0,004 0.006 0.008 0.01 
DU? 
^U-329291 
Fig. IV.1. Simplified expression for fractional release 


for the sphere. 


II 2 
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+ = Zei pe zt ZG) TAR LE)”. * (tV -2) 


Hence, D may be determined in a manner similar to that for the 
spherical case. For small periods of time, the correction factor 


3 pr ESE 1 =) is negligible; thus a straight-line plot may be ex- 
pected; NW, over longer time intervals the plot may be expected to 


show an increasing slope. 
F. Some Numerical Considerations 


Using the values of Shaked! for his sample $ 108 at 20409 C; ome 


obtains 


D=6.6x 10% cm/sec, 
t^-150 Ge9^ — (t- 2.85«10* sec), 


and G= 0.29 cm: 
Therefore, c. D = RI x ¡pe (CPs 154 «Jo ^) : 


The expected fractional release, with the use of Eq. (III-22), would be 


f- 5,7% = 52x10" 


Shaked estimated that the value of b he encountered during the 
anneal of this sample was 0.3 u/h. Because Shaked's sample was an- 
nealed in a Knudsen-type cell, this evaporation rate was considerably 
lower than would have occurred if the sample were completely exposed 
to vacuum. 


Thus, 


and, if the value of f calculated above is substituted into Eq. (IV-1), 


one has 





ae 


~3 p46 4-23 (bt) ,£ (01 = -3 
f Fe E 23 e IT es m m X lO 9 


with D representing a modified value of the diffusion constant 


D= tet ign cm/sec ) 


which is 61% less than that calculated by neglecting evaporation. 
If the sample had been completely exposed to vacuum at this 
temperature, the evaporation rate, b, estimated from thermodynamic 


dae would be about 150 u/h, from which 


6 OOOO 


For an anneal time of 2.25 X 10? seco and an actual D of 


DEDE E cm^/sec, it follows that 


"e, 
Vo m nm : 
/ ES p 
Using these data in Eq. (III-21), one predicts a fractional release of 
iuo Now, if the experiment had been performed in this way, 
and if the effects of evaporation were ignored, these values of f, t, 


and a used in the nonevaporation Eq. (III-22) would yield a value for 


D of 
ELT x pe G 


which is approximately 21,000 times as large as the actual diffusion 
eBeiii@tent for this temperature. 

These calculations illustrate the importance of either minimiz- 
ing evaporation experimentally, or of taking it into account analyti- 


cally whenever one obtains diffusion data at high temperatures. 
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V. CONCLUSIONS 


An equation has been developed that relates the fractional re- 
lease, f, of a tracer element in a geometrical body under conditions 
of evaporation of the body and diffusion of the tracer through the mov- 
ing boundary of the body, with these conditions occurring simultane - 
ously. This may be better expressed as the fractional release of the 
tracer element through a moving boundary. 

For small values of the quantity B, the release curve ís approxi- 
mately a straight line whose slope is slightly greater than that for the 
stationary-boundary case. As ß becomes larger in magnitude, the 
slope of the curve increases with increase in time and approaches a 
limiting value (see Fig. III.2) for the values of B and T considered. 

It can be seen that to completely ignore the evaporation of the 
body may lead to an incorrect value for the diffusion constant. The 
curve of f vs (1) 1/2 can appear to be a straight line when the bound- 
ary motion is not zero, but the slope of this line is not the same as 
that for the zero-boundary-motion line; hence, different values of D 
would be obtained unless a correction were applied to the moving- 
boundary case. 

Simplified expressions have been developed that, for small 


values of the variable Br 
determination of the diffusion constant, D. With these expressions 
one can predict a release curve whose slope, when f is plotted 
against the one-half power of time, will increase with time. The cor- 
rection made to the fractional-release equation for the stationary - 
boundary case is independent of D, and by subtracting this correction 
from the actual release, one should obtain a line of constant slope. 

It also should be possible to fit experimentally obtained 
fractional-release curves to the curves obtained in Sections II and III 
of this report, and thus to determine both the rate of boundary motion, 
b, and the diffusion constant, D. 


Experimental verification of the simplified and the approximate 


relations for the fractional release remains to be demonstrated. 
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Properties of the Error Functions and Related Functions 
Solution of the Slab Problem for a Moving Boundary 
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APPENDIX 
A. "Conventional'! Analysis of the Slab Problem 


Solutions to the heat-conduction equation for a slab -a<x<a, 


with a constant initial temperature To: and with T=0 at x=+a are 
given E 





Gigs E 2 A p l Cos ine DES (A-1) 


geo anti 


and 


TOT TEN ferfe ER ee a 


The-average temperature in the slab 1s given as 


O eo É (2n+l 
ty EL Zap ap [AGRA | a 


and 


Mr. Lo, E. Te "äer /) TE E 


Using an analogy to the heat-conduction solution, it will be shown that 
the concentration of a tracer element C(x, t) for any time t and any 
position x, where -a<x<a, the initial condition C(x, 0) 7 Cg, where 
C, is a constant, is 


0 
flee Ñ an«d)m^Dt N+ 
m 2 JC ZZ exp [- Las B (A-5) 
and 


LL See RSA ARRE L. 


E c 


Here, D is the diffusion coefficient for the particular element being 
studied. The corresponding fractional release, f, defined as the ratio 
of the decrease in average tracer concentration to the initial average 


concentration, is 


[- — Ca E) xi 6p E (A-7) 


Zw 


and 


f=2 


A.i. Tracer Concentration—Exponential Form 


le 





e E n. 
7 + 2EY verte A S (A-8) 





By using the method of separation of variables, it may be 


shown? that a solution for the diffusion equation 


2 C («) SC dë, (A-9) 


with C 6c o) ÓN (A-9a) 


ge Misato , (A-9b) 


may be expressed as C (x, t) - F (x)G. (t), (A. -40) 


where Feis A cosk x * B sink x ( A - 10a) 
n n n n n 


and G, (t) = exp EKFDE). (A -40b) 


Applying the boundary conditions (A-9b), yields B =0 and 


Ks a fi), where n=0, 1, 2,-°:. Applyimesthe initial condition 


(A-9a) yields a=Co, so Eq. (A-10) becomes 


E 


C, xt) - MOON exp (- Dt) 


(^ 2 
O h u > 3 | / | | l 


The coefficient A may be evaluated by setting t=0 and using the 


property of a Fourier sees that 


/ 
e OE EAS | JO 
An 77 f(x) COS = ax ; (A-12 


A 
thus, As Jl) cos Enel AX l where f(x) = 4. (A-13) 
=A 


This may be integrated directly to yield 


_,)" 
nico PL (A -14) 
" 77 (2n+1) 


Therefore, the substitution of Eq. (A-14) into Eq. (A-11) gives 


( (x,t) = Mr. 4, (-1)" exp J- fer JIDE eos rsrs ME 


m(2n+I) 


Any sum of solutions is also a solution, so the final result is 


C t)= E PR eL ener RE Las E (A-46) 


Bor long periods of time, the tracer concentration should become 


zero as t > æ., By noting that the exponent in each term approaches 


mt 


-o aS t-»o, it is clear that the tracer concentration will indeed ap- 


preach zero as a limit, 


A.2. Fractional Release —Exponential Form 


The fractional release, f, may be determined from 
Fus Q,— Q(t) (A-17) 
| or 


where Q(t) is the total number of tracer atoms per unit area contained 
is the total number of tracer 





mthe solidiat any time t, and Og = 2aCg 


atoms per unit area con Ue in the solid at time t=0: 


f. 246 - f (st) dx | 
M — 


=) OZ el ep Enero, Kun 


NC m 











= SI 2 nel) Fr Dt 
= 1-23 abp exp p cue m 


For to», the summation term approaches O as a limit; hence, the 


fractional release approaches 1 as a limit. For t=0, Eq. (A-18) be- 


Comes 
£= Je 8 = / (A-19) 
ES =: e. 
TE Flo (En+l 
VENE = 1 9 
The sum — + — +t kee = ————— is shown’ to be equal to 
4 32 52 n=0 (2n+1)2 


n 8 n^ 
—, so f=1-— X —=1-1=0 attime t=0, which is as expected. 


zoo 


Thus, it has been demonstrated how Eqs. (A-5) and (A-7) were 
obtained, and it has further been shown that the tracer concentration 
and the fractional release approach the proper limits for large periods 
of time, and at time t=0. 

It will now be shown, with the use of La Place transforms, that 
Eqs. (A-6) and (A-8) are the solutions of Eqs. (A-9) and (A-1417), re- 


spectively. í 


À.3. Tracer Concentration-—Error-Function Form 


Let oC (x,t) = D Set , (A-9) 
with C ( x, 0) 


and C (xat) UNIONS (A-9b) 


C. (A-9a) 


acto, =J (A-20) 


Now the La Place transform of C(x, t) -C(x, DS 
Eu o2 e 
C p)- f e-F Cogo - (A-214) 
O 


Making the transformations of Eqs. (A-9), (A-9b), and (A-20), one ob- 


by symmetry, 


tained 


eC (x,p)— Co = D LLLP), n) 


with D Es d. p) =o) NC S) 


and acto ) — O ^ (A-22b) 


02 


Defining m -p/D and by substituting in Eq. (A-22), one obtains 


2- Có p) 4 gC (iP) == — £e (A-24) 


Solving Eq. (A-24) for NES p), one has 





C tp) = Aef +Be + Ge (A-25) 


Upon applying Eq. (A-22a) and Eq. (A-22b) to Eq. (A-25), it follows 
that 


za a Colette) c 
C (x,P)= p (ef^-e7f^ = 


(A-26) 


Multiplying the first term top and bottom by podes one sees that 
St 
2 Co [536-9 en mmc C 
= — >=2 + =, (ASTE 


and upon expanding (1 E by the binomial theorem, one finds 


Hey] = a Zu on (ef - (e$). m 


(A-28) 
EI Hase 
A Ms 
Fis o 
By substituting (A-28) into Eq. (A-27), it follows that 
66 p)» - GE "T [ey CU] + 
HS P nzo P 
(A-29) 


A e [C2n+la-x] Ñ La Fe ye glen) +x] | 


N=0 


—- Co. Ce Seq 
d` Séi" 


sd 


Taking the inverse transforms, one obtains 


C£) G,- CE C" [enfe ERA erbo]. iaso 


For long periods of time, the tracer concentration should be- 
come zero as t--o, If t is allowed to approach infinity in Eq. (A-30), 


the erfc = > 1.0, and Eq. (A-30) becomes 


n € (ost) = Ce == ce ZCI") 


(A-34) 
- 6, - C,(e-2*2-2»--) 


which appears to be nonconvergent and equals C, or -C,, depending 


0 Dí 


upon whether the last term of the expansion is negative or positive, 
This is a periodic function, however; and if the function f(x) is defined 


as 
f(x) = 2 if 0 <x<r, (A-32a) 


and f(x) = 0 u -r SX S0, (A-32b) 


then the coefficients of the Fourier series for f(x) are 
O 77 
ro dX eL 2 dx =Z > (A-33a) 
TU 
6, f 2cos nx dy = I: Sin nx | cm (A-33b) 


m b> L[2 sin nx dx = 2000 nx] = E (1-c05 nr) ; (A-33€) 


and the Fourier series for f(x) becomes 


Ee por: 44.4 F (4, Cos nx +b, sin nx ) 


E 4» Sin Sx „au OX ru) = 
[us 2 = 


| 


A 





= / Mu, : L ony? (ua IX : ( A -34) 


On page 180 of Sokolonikoff and Roda Eni it may be seen that f(x) 

is well-defined on the closed interval -m €x €m except for the simple 
discontinuity at x=0, has a finite number of maxima and minima, 

and is defined for other values of x by the periodicity condition 
f(x+2w)=f(x); therefore, the Dirichlet conditions hold for f(x), and 


the Fourier series lor f(x) converges to 
L(x) = 5 [E T (A-35) 


for every value of x. 


In particular, at x=0, 


Pad PST eal, em 


Because the coefficient of the second term in Eq. (A-31) is shown to 
converge to the value 1, the tracer concentration will, in fact, ap- 
proach zero as a limit as t > œ. 


A.4. Fractional Release —Error-Function Form 


Calculating the fractional release from Eq. (A-17), one obtains 


f- 2aC,- f C Gt) dx 
Ae oue 


= Ia, (CES ero Lord (A-37) 


Sie N nt!) 
t Werke Bla I 


a 
LL (TB mula, Hu ld 
34 | ae ort Sehe + ¿Ne Ec dE 


aa 
which, with the use of Eq. (C-11), becomes 


+ = ES J [erf krg - ert Ëss 
- 7i ierfe Enea _ (er m" xx]! 


(A -38) 


= = (Zen [ierte > verfc "al 





=2(2 ot) erte o ët Tierken, Sr ert fg 1 


By letting n*1-2m, the last two terms may be combined, and by 
Heine, trom Eq. (C-7), that ierfc 0=7 1/2 one finds Eq. (A-38) to be 


f- 2(2 sf [77 h pae verte Eh | (A-39) 


For t= 0, the terms in Eq. (A-39) involving ierfe EE becornne cc 
pwesonly the first term E (PE)? remains. For t=0, this 
term and also the fractional release are zero, 

For large t, as t > œ, the terms involving the 


ierífc m. leric 0, where serien) S with these substitu- 
tions Eq. (A-39) becomes 


lm 1 = 5, (EE d AS | (A-40) 


By using the same analysis as that for the tracer concentration, it can 


be shown that the summation term converges to -1; thus Eq. (A-40) is 


lin £ = ng 


+ > 


aae 


which is indeterminate; therefore, it is necessary to refer to Eq. (A-37) 


to resolve the indeterminateness. REECH 


ff = LS tn [Berto Earle 
+26 A CHE 


erl SLR * ]dx (Aa 


D NEN 0 0 ll | 5 — ge EL. ONE 


j m * ` e, et a 
= 2 [RE dx = RECI) A 
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This previously was shown to be 1, so 


lim f=} l (A-42) 


Loin 


It has now been shown how Eqs. (A-6) and ( A-8) were obtained, and it 
has further been demonstrated that the tracer concentration and the 
fractional release approach the proper limits for large periods of 


time, and at time t=0, 


A.5. Correlation Between the Two Forms for Fractional Release 


Using the method of RC S one can show that Eqs. (A-7) and 


2 
(A-8) are identical for values of n < 14. Here, 


ee wn 








¿Dt | 
for ——< 1, Equation (A-8) becomes 
Ae 
a Cs ae (A-8a) 
mi 2 


sque 


Rewriting Eq. (A-7), one has 


p T DU T 2 
f- = EIE He —Z oh mee : mE || (A-43) 


The sum Wa ae ae — | + D EE (A -44) 
is 


M=- 009 


^ -xu2 
is, for small x, essentially f e du. Thus, 
= 00 


od 


-Xu* Ye 
Je Hie = (77%) (A-45) 
— Géi 
Br or 
where x for this case is m K is 1 or 4, such that 
4a 


MO uo E (va y^. re 


Now, if both sides of Eq. (A-46) are integrated over the limits 0 to x, 


the result is 


X + 25) oh (Le) = 2m” (A-47) 


which may be arranged 


= 


= x i 
BR, e té (X) Bs (A-48) 


Upon substituting Eq. (A-48) into Eq. (A-43), one finds, 


+ = [+ (reif - 2 As | 
-A fil-g -Ehl = 


m 


= |+£ ie Dt)” J^pt. Six | 
— Ttg a 30 Ke E 


Alr OI Dt) AS © (A-49) 








2 
Ja Zar Lim 
Te 
= -&27 
4 9 i 
The » — is shown’ as equal to er hence, 
sel m 6 
2 (Dt)” 
f= a (A-50) 
nZ Dt 
For values of > < 41, Eq.(A-8) reduces to Eq. (A-50); thus, for 
a 
mé Dt E 
these values of SCH Eqs. (A-7) and (A-8) are equivalent. 
a 


Solutions of the form of Eq. (A-7) converge slowly for small 
values of Dt/a* , but solutions of the form (A-8) are rapidly conver - 
gent for the same small values. ' This statement also applies to the 


equations for the tracer concentration, i,e., to Eqs. (A-5) and ( A-6). 


s 


B. ''Conventional" Analysis of the Sphere Problem 


Solutions to the heat-conduction equation for a sphere, O«r«a, 
with a constant initial tempe rature Tg with T-0 at r-a, and with T 


finite at r-0, are given as 


T(nt)- -s&te 7 EU ont, GH Sing UP 


n=] 


TOt) = -27 Že exp (- 288.) | =a 











e dla S&S Za LRT S 
D Tit) L lye LED erc ert | | 


Toter]. wm 


[See note following Eq. (B-6a).] 


The average temperature in the sphere is given as 


T, = HL Ea ep ei) Si 





and 


ly = Lë CT ZE erk ob | + 37, (GE 4E). (B-4) 


Using an analogy to the heat-conduction solution, one can show that 
the concentration of a tracer element C(r, t) for any time t and any 
position r, where 0< r<a and the initial condition C(r, 0) =C o 


with Co aon tant and with C(0, H) finite, i15 


-48- 





Cirt) =- E L erp ELA) sin BE (nas 
a 2 
R ea C,ZCU'exp (AE ; (B-5a) 
and 
0 SS Ja + ` 
= 0,382 rn EE, s 


C(o£) -G- E y LEE | (8-68) 


O 


[It will be noted that Eq. (B-6a) differs from (B-2a)--other than 
in the change of notation--by a factor of two in the summation term. 
It will subsequently be shown that Eq. (B-6a) is the proper expression 
for the tracer concentration at the center of the sphere.] 

The corresponding equations for the fractional release of the 


tracer element ame 


Pl REA) TN 


and 
M 6 (DF [r E verte ES | Fur 3( E) KEE 


D.1. Tracer Concentration —Exponential Form 


The radial-diffusion equation 


H Let, a ec 2 


-49= 


may be reduced to the linera-diffusion equation (A-9) by making the 
Substitution u(r, t)=rC(r, t) 


du Ic) c EN (UO, 
Ot aD (B-10) 
with the following boundary and initial conditions 
u(r, 0) = Tu I (B-10a) 
aa Tt Em (B-10b) 
and u(0,t) 20 E Beide) 


A solution for Eq. (B-10) may be obtained by using the method of 


separation of variable e 


Ur, t) = FA G (t) , (B-11) 
where F (r)=A cosk r+B_ sink r, (B-11a) 
n n n n n 
X ^ Dt 
and G (t) =ae (B-11b) 


Applying the boundary condition (B-10c) yields A= 0. The boundary 
nt 
a 
i. 072959 »*), and Eq. (B- 11) becomes 


condition (B-10b) gives K. = , whereas Eq. (B-10a) gives a =Co, 


U (r,t) = ps Sin ial exp (- k Dt) 
K BD. B, Sin DIE ex py n*m? DC 


( B - 12) 


The coefficient B may be evaluated by setting t- 0 and by using the 


property of a Fourier series that 


b, d fees em dx i (B-43) 


"IET 
Thus, 


A 
P. e c Csi? Br dr f where DECIES ( B-14) 


This may be integrated by parts to give 
PERO ed Se 
an ee Ie. 3 n = Ac ce DM (B-15 


The substitution of Eq. (B-15) into Eq. (B-12) results in 


^f 2 
U (rt) = au) CU Sg, rt) sin E : (B-16) 


Any sum of solutions is also a solution; so the final result, upon using 


Corn) zm de ES 


Me Zale FO exp |- P1 iat IEEE | sin ar yr (B E 
and 


C (oit)- pz Cnt)» 20, zc deep, - EMT (B-17a) 


where the relation 


lE 
x +0 x 


has been used. 

For long periods of time, the tracer concentration should be- 
come zero as t -»o. By noting that the exponent ingeachitermiin bor 
Eqs. (B-17) and (B-17a) approaches -œ as t = œ, itis clear that 


the tracer concentration will indeed approach zero as a limit. 


ze 


B.2. Fractional Release — Exponential Form 


The fractional release may be determined from 


i- (is c GU D (A-17) 


where Q(t) is the total number of tracer atoms contained in the solid 
at any time t, and Qo = = Ira” Co is the total number of tracer atoms 


contained in the solid pa time t= a Thus, 


f- fz ie — [te^ Cr 
u 1/3 Uc 


| - z& (- S26) z CU To p 9 dr md 


|- & E vh exp nem Pm 


Le 


For t+, the summation term approaches zero as a limit; hence, 
the fractional release approaches 41as alimit. For t-0, Eq.(B-18) 


becomes 


T |-& zx | (B-19) 


Ds a 
E Z 
The sum EE er ME ig shown oc re u 
4¢ 24 > 34 n=" F 6 
soe f=1- x e =4-4=0 attime t=0, which is as expected. 
T 


Thus, it has been demonstrated how Eq. (B-5) and ( B-7) were 
obtained. It has further been shown that the tracer concentration and 
the fractional release approach the proper limits for long periods of 
time and at time t=0, 

It will now be shown, by using La Place transforms, ' that 
Eqs. (B-6) and (B-8) are the solutions of Eqs. (B-9) and (A-17), re- 


spectively. 
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B 3. Tracer Concentration — Error-Function = orm 


We have 
GH, TR Z KH | (B-9) 
with Ga dE (B-9a) 
C(r, 0) 2 C, , (B-9b) 
and aC (af). o by symmetry. (B-20) 


The La Place transform of C(r, t) -C(r, p) is 
= yl » t 
O 


Making the transformations of Eqs. (B-9), (B-9a), and ( B-20), one ob- 


tains 


pC(np) -C, - D[ 22 cmp], 


(B-22) 
where C(a, p)=0, (B-22a) 
aC (Qo p) 
d Sn Se : B -22b 
Defining 
q^ - P/D (B-23) 


and substituting Eq. (B-23) in Eq. (B-22), one has 


akc (np) , 2 XNP) _ wi anne (B-24) 


OS 
Solving Eq. (B-24) for C(r, 5 one obtains 


pr A ed RE | gei 
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Upon applying Eqs. (B-22a) and (B-22b) to Eq. (B-25), it follows that 


1, .n-_Galet-et], Co. 
Se 


(B-26) 


Then, multiplying the first term of Eq. (B-26) top and bottom by 


- a e 
e d , one finds 


Cr ET el, E 2 s 29 





SE -1 


and expanding [1- e by the binomial theorem, one gets 


5 e = ia: e TA "ml Je E S (B-28) 
= < Ze 


Substituting Eq. (B-28) into Eq. (B-27), one sees that 
ch C Q C d Z 95 Ua -9 Lela C 


Upon taking the inverse transform, one finds 


eG) —C.— zZ lee - er | (B-30) 


and 





(C COT] < WE aian = au -£ > (B-30a) 


which is indeterminate. Applying L'Hospital's e to the second 


term to eliminate the indeterminateness, one sees that 


lim £C) _ lim f(x) 
Xa G(X) xoa q'(x ~ 


42 


Now, 


Cod)= Co- Im ASe Z [ere Corli- erf a | 
= C,~ ly, C (flle tT 
ier Ef iom 
ën in eo inte) 
-C.- Lalo EZ = EN | 


As an — check on Eq. (B-30b), it can be shown that the 





tracer concentration approaches zero as a limit as t >œ, 


The sum 

cx E 2 aD 

2080 ^w 1e e (A -44) 
N=- =] 


oo 


ISO sinaloa: TOS f MM 


du, 
o2 £ JL 
where fea = E) E ( A-45) 
— 00 


Equation (B-30b) may be written as 
X 2aC e 
4 (ot)=C, E ep Pe) Zal ane ^ (B-31) 
Combining Eqs. (A-44) and ( A-45), one finds that 


2 a 
Bee He A 
ze ÄIS) 2 ? 


(B-32) 


ED 


and by substituting Eq. (B-32) into Eq. (B-31), one obtains 


lo 
ie clos) €. rs. Ste Ae yi > | 


(ee (Tp SE 


Li UL E rs E (DE) Kell (B-33) 


C.-C. = 


Fractional-Release —Error-Function Form 


B.4. 


Calculating the fractional release from Eq. (A-17) 
E E 
a 2n4)a-r 
SKL cl Set, 
(B-34) 


—erfc& en 2 dr 
A[rÉ ferfe ST erf Geles dr. 


, one has 


n 


Now, Eq. (B-34) may be integrated by parts, by using Eq. (C-11) and 
by letting u=r and 


dv= erfc ce SEKR 


which yields 


f- & [a too ie ene. ip) erk Et 


zi d (Du orf ee T UIZ z Ze | = 


o 


= DC YR TS | na. =. n+1)a 
== GZ E Ler£c (pee + corta | (5:235) 
pais, na S. 
= 2( 23 Z Cerfc tbv. 2, ere DER 


By letting n+1=m, the dente tensis M5scOconmnib Mc tie 
¡Corfo terms will cancel excepirersam i^erfc O™term, Noting ther 
the ierfc Oe gl aZ orc OPER LEES TN 


comes 


f- Të [m "+22 (erfe (oa | psa. (B-36) 


For t=>0, the terms in Eq. (B-36) involving ierfc k/tt/2 become 
zero, Her t=0, the remaining térms are zene, thus, the tractional 


release is zero. 





Ihe “Ssubstitution,of BG.(G=7) 1s made for ieric in 


Eq. (B-36), the resulting expression becomes 
‘cs E (DE) Ve ) Im er A 2 exp - EZ) 
PHE 
-2 Z er cia | -S (ie) 


Horn largest zasıt >, > 0; thus, Eq. (B-32) may be substituted 


to Eda Hence 


x en fim aes tJ "ei - 22 B. erfe (oe 
+ roy] > -3 (2) | 
= lim | I2 Zlnderfe Bact 6(&)-3(&)| (B-38) 


— lim | 5(28)- PAG erte Dee | 


tr 


ERE 


Masis equivalent to œ= ©, which is indeterminate, Referring to 
Par B-34), one finds for large values of t that the difference be - 


tween the two complementary error functions may be approximated by 


T ee 7 HS lep - ST. Buen L sf | 
which, upon substitution in Eq. (B-34), gives 
lin f= Lin, X [Gi op [Er Ter 

= im S m x» [^ iz exp ma] -Zep 92 n'a SA: 


Applying Eq. (B-32), one finds 


im £ - 1 zn nor Er (Dt) — Hd zL-(rDt)* ldr 





/ / 
f irn RE ER] 
d 3 
= áj r = 22 = | 
THOS, m fads (B-39) 


It has now been shown how Eqs. (B-6) and (B-8) were obtained, and it 
has further been demonstrated that the tracer concentration and the 
fractional release approach the proper limits for large periods of 


mne and at time t=0. 
< 4, Eqs. (B-7) and 





2Dt 
Booth has shown that for values of E 5 


6 a 
(B-8) are identical. 


EO 


Solutions of the form (B-7) converge slowly for small values of 
Dt/a^, but solutions of the form (B-8) are rapidly convergent for the 
same small values, ' This statement also applies to the equations for 


the tracer concentrations, i.e., to Eqs. (B-5) and (B-6). 
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C. Properties of the Error Function and Related Functions 


The following definitions are made: 

2 f — 
Gp i. eme J e dz, (C-1) 

NIT o 
erk = ie (C -2) 
erf(-x) = -erf x, (C-3) 
9o 2 

erfc Ic E e’ dec (C -4) 


Nr X 


jerfc x = $ eric z dz, (C -5) 


X 


00 
iverfc x = f uod erie Z d2; (C-6) 
x 
10 
where i erfc x = erki. (C-6a) 
ierfc x = D - xXerfc x, (C-7) 
n 
and i“erfc x = : [erfc x - 2x ierfc x]. (C -8) 
The general recursion formula is 
m .n-2 .n-1 
DIU GNE = 1 erfc x - 2x i Erte. xe (C -9) 


Mota sides of Eq. (C=6) are differentiated, 


In d | En 
f = —— Nene d 
1 erice2 y ay 


Z 


d 
dx 
(C - 10) 


.n dz 
2 =] eric 72 5 
dx 


20s 


Integrating both sides of (C-10) and rearranging yields 


Du (can 
dz 


dx 





f iPerte z dx = - 


provided zZ is a constant and the constant of integration 1s neslected, 
then (C-11) is the indefinite integral of i'erfc z dx, where z isa 
function of x. 

The values of i"erfc( -X) will now be studied to obtain an inver- 
sion formula, so that the tabulated values of the positive argument 


may be used. 


Noting erf(-x) = -erf x (C-3) 

and ceri Fk =- Icen, (C -4) 
one obtains erfc(-x)-21- erf(-x)-2 1*4 erf x 

=4+(4-erfe x) (C -12) 


= 2 => eric x2 


1 x^ 


Thus, decas = e - x erfc x, (C -7) 
vn 


2 


e =) ae (=xjerfe (=x) 


and ierfc(-x) - 


e* 4 xerfc ( -x) 


SUE +x(2-erfc x) (C -13) 


A 
m 


2 
= 2x + E e ^ - xerfc x| 
Nm 


= 7X +lerie, x . 


= 


SO, pe Ses [erfc x - 2x ierfc xl, (C -8) 


A | 
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and deett (-x) -= L EAS eo (>) 
= > [(2 -erfcx) + 2x (2x+tierfc x)] 
j ; (C -44) 
= A = ASES x 2x Terie x) + =x 


1 


= > +x° - iverfc X 


By using the general recursion formula (C-9) and the relations just 


developed, it may be shown that succeeding terms are 


ee (-x) = : x + s - E D scc x. (C -45) 
4 4 1-2 4 4 .4 
eerie (=x) = — ger, ey serie x. C -16 
sl 16 4 12 \ 
5 1 1.9 deu 5 
and fc (-x) = — xt — x” + — i f Z C - 17 
n 1 erfc ( -x) m m £0 x i erfc x ( ) 


From the foregoing, a general inversion formula for i"erfc(-x) may 


be obtained: 


n o BH un zs | x74 
i erfe(-x) =(-1) i erfc x + —— T zen 
n! 2(n-2)! 
(C-18) 
n-4 n-m 
k E NET 2x 
16(n-4)! (4-8 -12-** 2m)(n-m)! 
where neme for Meo 15085501, d 
and H DS DU for nec 2 46, d x 


Dt 


D. Solution of the Slab Problem for a Moving Boundary 
D.1. Concentration of the Tracer 


An expression for the fractional release of a tracer element 


through the moving boundary of a slab will be obtained. Chambré S 


has derived a general solution for the tracer concentration for a slab 
-(a-bt) « x«a-bt and 0<t< = where b is the rate of movement 


of the boundary 


C (x,t) = zeen Ët op ESE Jo (D-4) 


ELP ep DE Gena) EIER 
O o cuo 


This equation is the solution of the diffusion equation 


9C(xt) _ D 2^ C (xt) , (D-2) 


Sueco the restricting conditions 


C (x,0) =P (x) tor O<X <4, (D-2a) 
OCI y fr ON (D-2b) 


and CGE NEN S (D-2c) 


For this problem, the initial condition (D-2a) shall be assumed to be 


a constant, Lo thus, 


C (x, O) = C. tar O«Ix|«4: (D-3) 


lf this sūübstitution and S(x t; £ F 2na) are made into Eg. (D-1) One 


has 


C(x t) = Siam |f exp -G£ 1a; + 


a 2 A 
serene ee o 
=A E 


This may be simplified somewhat by defining 





y =€ + 2na (D-5) 
and Z = uio ; (D-6) 
2 n Dt 


upon which Eq. (D-4) becomes 
CSR 
Cut) = Le Le | s exp (de 


d (DN) 


Ria ol a4) 


Ihe first term may now be integrated directly, and if exponentials are 


combined, one has 


C(xt) = $o ferf Leg, = erf [BeBe zl 
rs D _ 2G BO ]; E exp|- en rl 2 La nbe] (D-8) 
+exp|-* eg ab + Loy Dag: 
The exponentials may be simplified, 


T — EN 2 | 
bed! e = 7DE fx “ai #nbtyf 
-ohg (FO) eR Gareth = 


E 


= E Ir, -(x+2nbt) | — Hn bt (xenbt) } 





2 (D-9a) 
Se ly -er2nbt)] — nb (X +nbt); 
X F a 2 a | 
AO = dz [af bt 
qe lf V Aui 20H) e Gant) v x^-(X-Znbt) f 
(D-9b) 


=e [4 +(x -2nbt)] + 4nbt (x-nbt) \ 


Sr E Put x- -2nbt)] “ab (x-nbt) ; 


C6) - -£ er a - ert Ars] 


a(2ntl) E 
tyme SZ "exp|-2 nb(na-nbt- eiis (D-10) 


Lëscht Ptna-nbt+s)| m f. lico j^ 


te Së 


oimple substitutions similar to Eq. (D-6) allow the last two terms of 


Eq. (D-10) to be integrated: 


SEE 


fr 





d erf dios | 
+ £2 F¢i\"exp -ak (pa-nbt -x)| erf 28r- Grent) 


— erf &(£n-) -(X«2nbt | (D-11) 
2 ( Dt)" 


2 


vs ët "exp E nb(na- nbbex) erf 4 went! DI "n bt 


m sra 2nbt) ` 
e e 
Using Eq.(C-3) and rearranging terms, one has 


C (x,t) = E Jer Sim + erf Ztoey Zë, | 


zn exp LS ( na-nbt-x)] e r£ A 


er: | (D-12) 


ES Lia LR pande SREE 
T £n u)" SH | 


With the use of Eq. (C-4), i. e., erfe x=1-erfx, one finds 
A+X QX 
C(x t) = P = ferfe 2( „rerfc Zio | 
Cg c an L Ab (na-nkt X) lerfo en(a-E2-a.-X 
+ SE exp tere Re 


mE FOX 
-erfc co] | (D-13) 


z nb(na-nkt t de ríc pM 


For t=0, because erfc œ = 0, Eq. (D-13) becomes 













C. 
RES 














^£ 


C(x, 0) = Co . 


EOD 


For t-a/b, which corresponds to total evaporation, the tracer con- 


centration should become zero. Hence, 
C(x&.) - C, - & | erf GX, + ere SX, 
5 o R rs ZC) 
EE > 
NEE 


ZIN exp- njer Sie 2 Bay - erfc & SE zy | 


Since |x|<a-bt at t=a/b, x in Eq. (D-14) must be zero. Thus, 
(6, £)- T a G [Zerfc (30 ba. «e. ZE!) Dëst: 2a ST. 24 erc Ef If 
d Gë D erfc Gi N C E erte RA] = oti! 


with the use of Eq. (C- 12), 


- C, |I- -erfc z gl « o0, CECI [I-er Seel (D-15) 


oo 
It was previously shown in Eq. (A-41) that 2 » ane =1; similarly, 


00 n=0 

1t can be shown that 2 ES (-4)"= -1. With this substitution, 
n=1 

Eg. (Degen 


CE) C [Ier 29^] — C. []- ene 259)" | = Q- (D-16) 





Te 


a Ea. (D-13), b=0, which corresponds to no evaporation, then 


Gc.) = Co- $e o [erte Sx, t erfc Say | 
+52 ZN [er EH outen ` inan 


n- pax x 2n+l)a+X 
erk oS is 2Z(Dt)"z | 


+ 
If 2n-1=2m+4, Eq. (D-17) becomes 
2 A 
A en od 


AIR: = +i 5 
vere A preen + ork Bagless 


= 0 C Tp le Enea To SR A 





O 2 ño (DE A 
EN (ler D or «erc EmtDa.-x 
2m Ef A (Dt) 


tb Sie ees 


which may be recognized as Eq. (A-6). 


bee © Gactional Release of the Tracer 


The fractional release may be computed from Eq. (A-17), as 


follows: 
(a-bt) 
f- Q-A) _ 2al. -J C(xt)dx 
II E IB o‘ 
Qo 2ac, 


— 
—— 


SA 


y (0-bt) 


= Ls ue T GX t)ax ) 


"G-A 





using the symmetry relation, one obtains 


(a-tt) 


f= | | [ C69 dx (D-19) 
NS Ke E e S = 
ato y 


With the subst: Mob of Eq (D- mu. lor coU x. t) Vie follows tae 


bt) 
E (= a fe, d ‘ate Co ferke Ae RN Id 








e [Eze Vex p[-Ba(na-nbt- der: au bo as X 
` E (D-20) 
GI 


-z JL P zc I) exp nb (na-nbt+x) erre. m Hl o 
- erc AG Eg ëch dx. 


With the use of Eq. (C- 11) in the second integral, Eq. (D-20) 


becomes 


c Zee eje nb orf ae -X 


- ero a ] dy 


(&- bt) 


5 Z C exp E (nya le - nox.) [erf Ent Rna Eon AX 


a T = 





ea 


Po =bt 


2 
F. "2 l aro da ero TA m 
Zot tol en cea 


@-bt) 
= aen als ae a be erf SE 2 kt 
O 


The remaining integrals in the expression (D-24) for í may be 


evaluated by integrating by parts, as follows, Let 


= exp (22x) ad — dy-erfc en, L ES 


du Ce and \/= (DY ierfe net a 


The result may be written as 


eo A o TO Ver meer 2 


To integrate the second integral by parts, let 


2n(a-bt)s d L N | 
y= expt zex) and dy- erc tr : then 


du s Mb ey 1d and y=-2(D) Pce rfe GE RTL 


E 


The result is 
E b^ XC /) "expt A BAR ao: E in ee gtk (D-23) 


Upon applying the limits to Eqs. (D-22) and (D-23) and inserting the re- 
sults into Eq. (D-21), one obtains 


f= BE bt)" ect Aë Zë erte Sëch, 


/ Lo 02 yt! qx ml 
De ZU Tev Lë) Ze I) pey] h p pe H) rf Se 
—("erfc EE — €Xp Be) ree Desde (D-24) 
41 Uere Aosta en xpF- Oh Alert EE E 


: i g ra. 
Uert: Se TOP, ln SE RE oer 


After eliminating equal and opposite terms and combining exponen- 


tials, it follows that 


f- «(ey lierfe ee, ¿enfe SKF, | 


fr) Ze thon. HF zi" pato "eroe ail a 


-apta "erk Er atta t exp EJ "rfe nella HU jo 


ROC) TR ADO) V2, 





Lie 


i IYa-bt)-a 
- exp Dë be) | (nae ER C 


MT. ; -6C 
= Ht + (Bt) ent ap = ierk ERA 


SIS Ze DA (oT: n LAN, yest L C 
E ‘marten: zc ye 


«(E J. EC (EG br! exp |- oves in hr Gg SS 


—( Morfe Anl a-bi)-a. ; 
2 (D6) 7> 
At t=0, each term of Eq. (D-25) becomes zero; thus, the fractional 


release is zero, Attime, t-a/b, which corresponds to complete 


evaporation of the solid, the fractional release should be 1, as shown: 
2: Q . De 
E + lierfe 375 (éen ` certo faja | 
nem[ - rnb Y DaY2 im BINA HN 
d y: zl Loës y 1 erte Zr Bay 
, A : a. a Er 
— (Moyfc ZC la a (orl Spe (Merfe Aa | 


If, in Eq. (D-25), b is set equal to zero, which corresponds to ro 


(052) 


evaporation, then 


f= (DN Tuerto SE Lertc Ge EF: |+ HE 35; ZI" enc De 


(D-27) 


c na ntl)a. ` 
Le rfe 7% (Do 4 ceríc DDR mo sch 3 


E 
by letting n-4=p and n+1=m in Eq. (D223), enesonitams: 


(PY lierte o -ierk gya + RCW ere Bh. 


E | eo +2. 
Lë ierk d, "SCH certo ty. | 
(D-28) 


= LS ierfe O 1 Sëtl verte n 


Lae + ZN erh, 


Noting that ierfc osr 2 and combining the m and p terms, one 
finds 
T > 
ah 
f-2& T ^4 2 CI) cere d (D-29) 


which may be recognized as Eq. ( A-8). 


D.3. Conversion to Dimensionless Form 


If one defines 


T Z — (D-30) 

ae 
and B = = | (D-34) 
so that E E La E EE. a «gt «! À (D-32) 


where QgT-0Ü corresponds to t=0, and Br=1 corresponds to com- 


plete evaporation, then Eq. (D-13) becomes 





rfc AR aro 5% x e 4 





az 














i - en 2n s: p» d 4] 


2v" 


Equation (D-25) becomes 


+=pT+ Tr lief Eur P ze —_ ) 


n n (- Y ten T SE L Renee - 1 


P) AS 


=a mo = Rn- cl 
e era Ss (D-34) 
^ Le 


leue Yn, mo 
rt EN (¿nac ^J exol- ne (nix! 49 uM E 
— merlot SE Tei |: 


D+] 





iJ 


and Eq. (D-29) becomes 


de E = T Kì 
f= E E ¡7 J- TeC o : (Dea 
| E 
D.4. Approximate Solution for the Fractional Release 


If Eq. (D-34) is expanded to the first few terms in n, then 


f=B7 po SC orto fox -ierk £ =) 


- t^s. We (96 zech! Ne x nord, E TX 


OS 


EN YT aaa adla SER inen 2382 
+ UZ Cape’ exp [2-40 inerte 288, et ZE en 
VU ec lap alee erte SAE - erk LE) 
Aë Jy" (eec Je exp E Gefeiert Za t- merk Z2 Se 
ZG Ger exp a] (Perf GER - "eric Sr) 
en. 


Now, if the restrictions d < 0,01 and 67<0.5 are made, then 
a 
Eq. (D-36) may be closely approximated by 


f=AT+T Perk E er ABE cre Xt . (0-3 


Inversion formulae for i''erfc (-x) have been developed in Appendix C. 


If these are applied to Eq. (D-37) and expanded in terms 2% itcan 
be seen, if sufficient are taken, that terms of ps y cancel 
and only the (erte Er 





terms are lckt Thereupointhe approxi- 


mate expression for the fractional release, f, becomes 
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E. Solution of the Sphere Problem for a Moving Boundary 
Pear Concentration of the Tracer 


In this section an expression for the fractional release of a 
macer element through the moving boundary of a sphere is obtained: 
Chambré : has derived a general solution for the tracer concentration 
fora sphere O« r«a-bt and O«t« a/b, where b is the rate of 


movement of the boundary 
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This equation is the solution to the diffusion equation 
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If this substitution and S(r, t; £ + 2na) are made into Eq. (E-1), then 
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Integrating the second term of Eq. ( E-5) with the aid ot Eq. (6-11 
and also substituting the results of Eqs. (E-6a) and(E-6b), one finds 
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The remaining ıntegrals may be solved as before by integrating by 
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Upon integration of the last two terms of (E-8), one has 
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By using Eq. (C-12), one obtains 
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Using Eq. (C-13), one has 
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The substitution of Eqs. (E-10) and (E-11) give 
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For t=0, since the erfc o = 0, Eq. (E-13) becomes C(r, 0) = Co: 
For t=a/b, which corresponds to total evaporation, the tracer con- 
centration becomes zero. (The proof will not be given here.) 
If, in Eq. (E-13), b=0, which corresponds to no evaporation, 
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which may be recognized as Eq. (B-6). 
E.2. Fractional Release of the Tracer 
The fractional release may be computed from Eq. (A-17), 
f = —, (A-17) 
G 


where QO(t) is the total number of tracer atoms contained in the solid 
at any time t, and Qo - i EEN is the total number of tracer atoms 


contained in the solid at time t=0. 
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The first integral in Eq. ( E-17) may be handled directly, whereas the 
second and last integrals may be handled by integrating by parts twice, 


as follows with the use of Eq. (C-41). 
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The remaining integrals in the expression (E-18) for the fractional re- 
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differs from Eq. (E-19) only in the power of the ierfc x term. If m 
were replaced by m+1, the result [Eq. (E-19)] could be used as it is. 
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Upon applying limits to Eqs. (E-19) and (E-20), the resulting terms are 
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as the solution in full form for the fractional release of the tracer 


element {rom a spherical body. 
At t=0, each term of Eq. (E-214) becomes zero; thus the frac- 
tional release is zero, Attime t-a/b, which corresponds to com- 


plete evaporation of the solid, the fractional release should be 1, as 


follows: 
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By letting n+4=m and n-4=p, one obtains 
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Because ierfc 0= -1/2 and i^erfc 0- 1/4, Eq. (E-24) becomes 
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which may be recognized as Eq. (B-8). 
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E.3. Conversion to Dimensionless Form 


If the definitions (D-30) and (D-31) are now utilized, along with 


Bart D-32), then Eq.(EH-13) becomes 
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Equation (E-21) becomes 
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and Eq. (E-25) becomes 
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E.4. Approximate Solution for the Fractional Release 


Now, if the restrictions => <0.01 and 67<0.5 are made, then 
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Eq. (E-27) may be closely approximated by 
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Application of the inversion formulae for il erfc( -X)--developed in 


Appendix C--to Eq. (E-30), upon expansion to a sufficient number of 
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terms, reveals that terms in ag cancel each other. There- 


upon the approximate expression for the fractional release will reduce 
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In general, the last three summations of Eq. (E-31) may be neglected 


in comparison with the remaining terms 
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NOMENCLATURE 


Characteristic dimension at initial conditions; half thickness for 
a slab; radius for a sphere (cm) 

Rate of boundary movement (cm/sec) 

Concentration of the tracer element (atoms/cm”) 

> 


Diffusion coefficient for the tracer atoms (cm* /sec) 


Initial concentration of the tracer (atoms/cm 


Fractional release of the tracer element 

Number of tracer atoms 

Initial value of Q 

Space coordinate in spherical geometry (cm) 

Time from start of evaporation/diffusion (sec) 

Space coordinate in slab geometry (cm) 

Dimensionless evaporation/diffusion coefficient Beg (II-8)] 
Integration variable (cm) 


Dimensionless time [Eq. (II-7)] 


10. 


EE 


REFERENCES 


Hagai Shaked, Diffusion of Xenon in Uranium Monocarbide, 
UCRL-10462, Nov. 1962. 

G. Long, D. Davies, and J. R. Findlay, Diffusion of Fission 
Products in Uranium Dioxide and Uranium Monocarbide, 
TDCI Oct 1960. 

J. Belle, Uranium Dioxide: Properties and Nuclear Applica- 
tions (U.S, Atomic Energy Commission, Washington, D. GC. 
July T1961): 

L. B. Loeb, The Kinetic Theory of Gases (Dover Publications, 
leme New ork, 4961). 

P. L. Chambré, On the Diffusion of a Radioactive Tracer 
Element Through a Moving Boundary (published). 

A. H. Booth and G. T. Rymer, Determination of Diffusion Con- 
stant of Fission Xenon in UO, Crystals and Sintered Compacts, 
AECL-692, Aug. 1958. 

HS Carslaw and J. CG Jaeger, Conduction of Heat in Solids 
(Oxford University Press, London, 1959). 

RAM Barrer, Diffusion in and Through Solids (University 
Press, Cambridge, England, 1941). 

J. S. Sokolnikoff and R. M. Redheffer, Mathematics of Physics 
and Modern Engineering (McGraw-Hill Book Company, Inc., 
New York, 1958). 

E. S. Smith, M. Salkover, and H, K. Justice, Unified Calculus 
(John Wiley and Sons, Inc., New York, 1947). 





This report was prepared as an account of Government 
sponsored work. Neither the United States, nor the Com- 


mission, nor any person acting on behalf of the Commission: 


A. Makes any warranty or representation, expressed or 
implied, with respect to the accuracy, completeness, 
or usefulness of the information contained in this 
report, or that the use of any information, appa- 
ratus, method, or process disclosed in this report 
may not infringe privately owned rights; or 


B. Assumes any liabilities with respect to the use of, 
or for damages resulting from the use of any infor- 


mation, apparatus, method, or process disclosed in 
this report. 


As used in the above, "person acting on behalf of the 
Commission" includes any employee or contractor of the Com- 
mission, or employee of such contractor, to the extent that 
such employee or contractor of the Commission, or employee 
of such contractor prepares, disseminates, or provides access 
to, any information pursuant to his employment or contract 
with the Commission, or his employment with such contractor. 























iit itii j 





